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^1 HEX Optical spectra covering the evolution from about 6 days before to 

about 5 weeks after maximum light and the ROTSE-IIIb unfiltered light curve 
^ of the "Branch-normal" Type la Supernova SN 2005hj are presented. The host 

' galaxy shows H II region lines at redshift of 2; = 0.0574, which puts the peak 



C/3 



unfiltered absolute magnitude at a somewhat over-luminous —19.6. The spectra 
show weak and narrow Si II lines, and for a period of at least 10 days beginning 
around maximum light these profiles do not change in width or depth and they 



O 

j-j ■ indicate a constant expansion velocity of ~ 10,600 km s ^. Our observations 

indicate that Si is confined to a relatively narrow velocity region. We analyzed the 
observations based on detailed radiation dynamical models in the literature. The 
models considered include delayed detonations, deflagrations, pulsating delayed 
' detonations, and tamped detonation/merger scenarios. Whereas the first two 



classes of models have been used to explain the majority of SNe la, they do 



• not predict a long velocity plateau in the Si II minimum with an unvarying 

10 ' line profile. Pulsating delayed detonations and merger scenarios form shell-like 

^ ■ density structures with properties mostly related to the mass of the shell, Mgheii, 

O ■ and we discuss how these models may explain the observed Si II line evolution; 

however, these models are based on spherical calculations and other possibilities 
^ ■ may exist. SN 2005hj is consistent with respect to the onset, duration, and 

■ velocity of the plateau, the peak luminosity and, within the uncertainties, with 

the intrinsic colors for models with Mgheii = 0.2Mq. Our analysis suggests a 
distinct class of events hidden within the Branch-normal SNe la. If the predicted 
relations between observables are confirmed, they may provide a way to separate 
these two groups. We discuss the implications of two distinct progenitor classes 
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on cosmological studies employing SNe la, including possible differences in the 
peak luminosity to light curve width relation. 

Subject headings: Supernovae, SN 2005hj, cosmology 



1. Introduction 



Type la supe rnovae (SNe la) are t hought to be thermonuclear disruptions of white 
dwarf (WD) stars (iHoyle fc FowlerlllQGOl ). but the details remain uncertain. One possibility 
for the progenitors is the single degenerate model in which main sequence stars or post 
main sequence red giants transfer mass to a WD through Roche Lobe overflow or a common 
envelope and the WD grows close to the Chandrasekhar mass Mch- Due to compression, 
the thermonuclear runaway starts near the center leading to the explosion of the WD and 
a rapidly expanding envelope with a mass close to Mch- A second possibility is the double 
degenerate model in which a pair of WDs merge and lead to an explosion. In most such 
cases, the resulting mass of the rapidly expanding envelope will be different from Mch- 
Residual material from these mergers surrounding the explosions will get swept up by the 
ejecta forming dense, shell-like structures. Although the explosion of single WDs seems to 
be favored for the majority of objects, we may expect mergers to contribute to the SNe la 
population (see [|2]). 

The possibility of different progenitor channels, the population of which may vary with 
redshift, may pose a challenge for the use of SNe la in cosmological studies that rely on 
a single parameterization, such as a light curve width to peak luminosity relation, LWR, 
to reduce the intrinsic scatter in the peak magnitudes and render them standard candles 
( lPhillipslll993l : iPerlmutter et al.l 119971 ). To first order, the LWR relation can be understood 



as a result of differ ent amounts of ^^Ni produced during the explosion (IHoflich et al.lll996 



Umeda et al.lll999l ). There may be some spread and an offset in LWR introduced by one of 
the channels if the masses of the envelope differ from Mch, and/or the density structures 
differ. This can lead to a systematic shift of LWR with redshift if the evolutionary time 
scales of the progenitor systems differ. Even if the different progenitor scenarios obey the 
same LWR, differences in the color could introduce systematic errors in cosmological studies 
because SNe la are known to suffer to some degree from reddening in their respective host 
galaxies which has to be taken into account. To correct for this, the maximum light color 
excess (usually Eb-v) and an average reddening law are used to determine the amount of 
absorption. SNe la that are intrinsically redder as compared to the average local sample will 
thus be over-corrected in this fashion to a higher luminosity. 
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Similarly to the two distinct progenitor channels, qualitative variations in the explosion 
physics may lead to various classes of SNe la even within the single degenerate scenarios. 
Standard explosion models include delayed detonations (DD) and deflagrations. In these sce- 
narios, burning during the deflagration phase leads to an unbound WD. In DD models, the 
deflagration turns into a detonation in an expanding envelope. Because the density structure 
of the WD declines monotonically with radius, the resulting density structure in the expand- 
ing envelope also smoothly declines with mass and radius. A variat i on of DD models ar e 
the pulsating delayed detonation models (PDD; iKhokhlov et al.lll993l : iHoflich et al.lll995al ). 
In these models, the total energy production during the deflagration phase is, by construc- 
tion, lower and insufficient to unbind the WD. This results in large amplitude pulsations. 
Because the fall-back time increases with distance, the inner regions contract and leave a 
shell behind at larger distances. Due to infall driven compression, a detonation is triggered, 
the material behind the burning front is accelerated, and this expanding material runs into 
the low-velocity material left at larger distances. Similar to the merger scenario, a shell-like 
structure is formed with very similar light curve a nd spectroscopic properties, but with a 
total mass close to Mch ( IHoflich fc Khokhlovlll996l ). 



These two groups, consisting of DD and deflagration models such as W7 (INomoto et al. 



19841 ). which lack shells, and the models with shells (mergers and PDDs), can be dif- 
ferentiated by their p r edictions for the photosph eric evolution and maximum light colors 
( IKhokhlov et al.lll993l : IHoflich fc Khokhlovlll996l ). For the former group, the photospheric 
velocities, Vph, smoothly decline with time and the models show a blue color at maximum 
light, B — V ~ 0*"; in the latter group, Vph shows a plateau in the evolution as the photosphere 
recedes through the shell. These models are intrinsically redder and slightly over-luniinous 
because of the lower expansion rate in the inner region. As shown in IKhokhlov et al.l (119931 ) 
the color, length and velocity of the plateau are correlated with the mass of the shell, and 
this potentially allows the two groups to be distinguished even for similar brightnesses. 



In deed, there is a growing sample of SNe la s 



1990N; 
1999ee: 



Hamuy et al.l 



eibundgut et al.lll99ll : Mueller fc Hoflich 



2002 



lowing photospheric yelocity plateaus (e.g . 



19941 : 1991 T, 1999aa: iG aravin i et al.ll2004 



2000cx; iLi et all 1200 lal : see also iBenetti et al.ii2005l ). Many of 
these SNe la have been reported as having a red color B — V at maximum, but this is typically 
attributed to reddening along the line of sight. Alternatively, this sample may suggest the 
contribution of events with shell-like density structures in the observed population. These 
events may be understood in terms of mergers or PDDs; however, the inhomogeneities and 
incompleteness of individual data sets in the literature preclude definite conclusions. 

To address this problem and others, we started the Texas Supernovae Search (TSS; 
Quimby et al. in prep.) with the goal of providing a homogeneous set of quality data for 



-4- 



several supernovae beginning well before maximum light. In this paper, we present our 
observations of SN 2005hj and analysis of the data. In [|2]we describe the discovery and give 
the details for both the photometric and spectroscopic follow-up. In ^ we discuss generic 
properties of explosion models and suggest a secondary parameter to separate models with 
and without shells, and analyze the peculiarities of SN 2005hj. Conclusions and discussion 
are presented in §11 



2. Observations 



SN 2005hj was discovered on October 26.13 UT in the field of Abell 194 as part of the 



TSS. The TSS uses the wide field (1?85 x 1?85) 0.45m ROTSE-IIIb telescope flAkerlof et al. 



20031 ) at the McDonald Observatory in Texas to scan nearby galaxy clusters nightly for tran- 
sients with a modified version of the PSF-matched image subtraction code from the Super- 
nova Cosmology Project. SN 2005hj was found at an unfiltered magnitude (calibrated against 
the USNO-Bl.O R2) of Cr = 17.4 and is located at a = 0l'^26"^48!27, S = -0 1°14l6^^8. The 
foreground reddening at this location is Eb-v = 0.039™ ( ISchlegel et al.lll998l ). Examination 



of ROTSE-IIIb images from Oct. 20 and Oct. 22 shows the SN was detected prior to discov- 
ery, but not significantly well to pass the search pipeline's automatic cuts. Figured] shows the 
ROTSE-IIIb light curve for SN 2005hj through 40 days after maximum light. To construct 
the light curve, we co-added images taken on a given night (usually 6) excluding any frames 
of significantly lower quahty due to passing clouds or wind sheer, and then subtracted the 
reference image convolved to the same PSF. Magnitudes were determined by fitting the local 
PSF (derived from the co-added nightly images) to the location of the SN on the sub tracted 
frame usin g custom softwar e and the DAOPHOT PSF-fitting routines (IStetsonlll987l ported 
to IDL bv lLandsmanlll989[ l. 



The unfiltered CCD response of ROTSE-IIIb has an approximate full width of ~ 4000 A 
centered in the i?-band around 6000 A. Because we do have some sensitivity in the blue and 
since the B — V colors of SN e la typically grow ^ 1.0 mag redder in the 30 days after 
maximum (jPhillips et al.lll999l : iKrisciunas et al.l 120031 ). there is a blue deficit at later times 
that causes our unfiltered magnitudes to decline more rapidly than the true i?-band fading. 
Note that V — R colors of SNe la are close to zero at maximum light. We therefore limit the 
light curve fitting to data taken before 10 days after maximum (determined through several 
iterations of t he fit), during whi ch the color evolution is minimal. The best fit i?-band 
template from iKnop et al.l (120031 ) is also shown in Figure [TJ The date of maximum light 
determined from the fit is Nov. 1.6 with a formal error of 0.7 days (note the template ph ases 
are relative to the i?-band maximum). The best fit stretch factor (jPerlmutter et al.lll997l ) for 
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the light curve width is s = 1.2 ± 0.1. The prehminary measurement of the observed B — V 
color at V maximum from the Carnegie Supernova Project is 0.07"^ ± 0.05 after removal of 
the host light but before any extinction or fc-corrections are applied (M. M. Phillips, private 
communication) . 

Near real-time photometric analysis combined with target of opportunity (ToO) time on 
the neighboring 9.2m Hobby-Eberly Telescope (HET) allowed us to obtain optical spectra 
just 4 hours after the discovery images were taken and every few days over the next 6 weeks. 
These observations are detailed in Table O The instrumental response is such that very 
little second order light is expected blue of 8900 A even with the GG385 blocking filter. 
The data were reduced in the optimal manner using and custom IDL scripts. The 

wavelength scale was calibrated against Cd and Ne lamps and its accuracy was verified by 
comparing night sky lines to the spectral atlas of iHanuschikl (120031 ). Because the HET pupil 
size varies for different tracks, absolute flux c a librat ion ca nnot reliably be achi e ved; h owever, 
we used the standard stars of iMassey et al.l (119881 ) and iMassey fc Gronwalll (Il990l ). which 
were observed using the same setups, to achieve relative spectrophotometric calibration and 
to remove strong telluric features. 

The redshift of the host galaxy was derived from narrow emission lines around 7000 A 
(observed), which we attribute to H-a, [N II], and [S II] in the host galaxy. We combined all 
the spectra and simultaneously fit these lines with Gaussians to determine the line centers. 
The line redshifts are best fit by 2; = 0.0574±0.0002, and we adopt this value for the SN. This 
gives SN 2005hj an absolute peak magnitude of —19.6 in our unfiltered band pass (assuming 
Ho=71 km s" -^ Mpc~^, VLm = 0.3, an d Q\ = 0.7), and places the host well behind Abell 194 
{z = 0.0180; IStruble fc Roodlll999l ). The brightness and broad light curve shape suggest 
that SN 2005hj is a slightly over-luminous SN la. 

The unfiltered ROTSE-IIIb reference image shows that the host for SN 2005hj is rela- 
tively bright {Cr = 17.8) and compact, and is therefore likely a significant contaminant to 
our spectra. Thus, we have to subtract the galaxy contribution (see Fig. [2]). Lacking an 
observed spectrum for the host galaxy excluding the SN light, we constrained the galaxy 
SED using archival Sloan Digital Sky Survey (SDSS) ugriz observations and obtained a 
template galaxy spectrum (N. Drory 2005, private communication). The relative amounts 
of SN and galaxy light in the spectral apertures will vary not only with the changing SN 
brightness, but also with the seeing, slit width and positioning. Also plotted in figure [2] 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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is a spectrum of SN 1999aa (blue cur ve) constructed via a linear interpolation of the —7 
day and —3 day spectra presented by iGaravini et al.l (120041 ) . Noting the similarity of the 
spectral features of SN 1999aa and SN 2005hj, we assume that we can model our observed 
spectra as a linear combination of our galaxy template and the SN 1999aa spectra inter- 
polated to the same phase as the SN 2005hj observations. We perform a least squares fit 
to determine the relative contributions of each component. The red line in figure [2] shows 
the derived contribution of galaxy light in the —6 day spectrum. Aside from a few small 
differences (most noticeably in the Mg II A4481 triplet), some of which may be explained 
by calibration errors, the combined SN 1999aa + host spectrum (purple curve) is a good 
fit. The over all fit is improved if we interpolate the SN 1999aa spectra to —5 days instead 
of —6, especially in the 5400 A to 6500 A range, which could imply a ~1 day error in the 
date of maximum light or different time scales for the spectral evolutions of the two SNe. 
We repeated this process for all the SN 2005hj spectra, each time using the same galaxy 
template and the SN1999aa spectra (interpolated to the appropriate phase) as reference to 
determine the relative amount of galaxy light. In general, the galaxy template added to the 
SN 1999aa spectra does an excellent job of reproducing the observed SN 2005hj spectra. 
The galaxy light typically dominates the fiux red of 7000 A. Figure [3] shows the spectral 
evolution of SN 2005hj recorded by the HET between days —6 and +34 with the derived 
galaxy contribution subtracted. 



2.1. Spectral Characteristics of SN 2005hj 

Overall, SN 2005hj shows spectra with lines dominated by intermediate mass and iron 
group elements as is typical for SNe la. While the lines show normal expansion velocities, 
the absorption components are more narrow and, for the early phases, weaker than typically 
observed, as exemplified by the Si II A6355 line (see Fig. Hj). SN 2005hj also shows an 
atj^ical velocity evolution of these features over time. 

Line minima are useful diagnostic indicators of the ejecta structure as they give the 
abundances and velocities of the material. The actual measurement of the velocity at the 
minimum of the line profile is complicated by the presence of the continuum, other blended 
lines, and some uncertainty in the true line profile shape. Detailed modeling is required to 
accurately sort out all the components and how they relate to the photospheric layer to reveal 
the velocity distribution of the Such models have shown that the absorption minima 



^In general, Si II lines form above the photosphere and velocities measured from such absorption minima 
can be 1000 to 2000 km s~^ larger than those measured from weak lines. However, for shell models the steep 
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appro ximate the photospheric expansion veloc ities to within about 1000 km s ^ at maximum 



hght (IHoflich et al.lll995al : iMarion et al.ll2006l ). Thus, simple line fitting can lead to a rough 
description of the ejecta velocities, and allows a useful comparison to discriminate between 
different models. At late times the photosphere will recede below the Si rich layer and so 
the velocities derived from the Si II A6355 line will become increasingly discrepant with the 
photospheric velocity. For deflagration and cla ssical detonation models, t his departure will 
begin to set in 1-2 weeks after maximum light (IHoflich &: Khokhlovlll996l ). 



The strength of the Si II A6355 line and its persistence from at least 2 weeks before to 
4 weeks after maximum light make it a valuable tool for probing the ejecta. Its evolution 
with time is shown in figure [5] for the case of SN 2005hj. To determine the velocity, we 
smooth the spectra by Fourier Transform filtering, divide by the estimated continuum, and 
then select the lowest point using spline interpolation over a selected range. The continuum 
is represented by a medium order (6th - 7th) polynomial fit to regions of the spectra that 
are not strongly affected by lines. To smooth the spectra, we use a Fourier Transform to 
convert the data into a power spectrum, and then multiply this by a filter to remove high 
frequency variations. We then apply a reverse FT to the filtered power spectrum to recover 
the smoothed spectrum. The filter has the functional form 



f 1 for - z/c < 1/ < z/c /-j^N 

\ exp[— (z/ — z/c)^/2(T^] otherwise 

The filter cutoff frequency, z/c, and attenuation scale, a, were determined as follows: 1) the 
spectra were converted into a power spectrum, -P(z^), via Fourier Transform; 2) the slope 
of log(P) is fit over the noise dominated high frequencies and interpolated through the low 
frequencies to determine the noise spectrum; 3) z/^ is taken as the frequency at which log(P) 
drops to within three times the dispersion about the noise spectrum; 4) a is chosen such 
that the slope of \og[P{Pc + 2(t)] is twice the noise spectrum slope (i.e.. z/ = z/c + 2(T is the 
frequency above which noise is clearly the dominate component). For this analysis, only the 
spectral bins with signal to noise above 25 were considered (note the peak throughput for 
HET/LRS is near the Si II A6355 line). For consistency, we adopt a single filter for all our 
analysis, choosing the results from our nosiest data, Vc = 0.0066 and a = 0.0053 A^^, 
which removes noise in the data but also some real information related to "sharp" features in 
the spectra such as the narrow core to the Si II A6355 absorption in the day +10 spectrum. 

Using the relativistic Doppler formula and the /-weighted Si II A6355 rest velocity in 
the host galaxy frame, we convert the wavelengths of the line profile minima into expansion 



$z/ 



density gradients cause even strong lines to from very close (in radius) to the actual photosphere. 
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velocities. For each spectrum we conducted 250,000 Monte Carlo simulations in which nor- 
mally distributed noise based on the statistical flux errors was added to the data and the 
FT smoothed minimum was found. The peak of the distribution and the interval containing 
68% of the simulation results were used to calculate the velocity of the minimum and its 
error, respectively. We also measured the relative shift in the H II region lines over all epochs 
and found the scatter to be 80 km s~^, which we add in quadrature to the individual errors. 
The results are given in table [2] and plotted in flgure [61 We find that the data points are 
at 10,600 ± 150 km s~^ between maximum light and +18 days, somewhat faster prior to 
maximum, and significantly slower on day +25. By day +34, the Si II A6355 absorption has 
all but completely disappeared. 

From maximum light through day +10, the Si II A6355 line profile shows little change in 
both depth and width in addition to maintaining a constant absorption minimum velocity. Of 
specific relevance is the blue wing of the absorption profile; this section of the line is formed 
by the material at the greatest distance from the photosphere and at the highest velocities, 
and as such it should be the first to vanish as the photosphere recedes. The consistency of this 
blue wing from maximum light through day +10 suggests the photosphere falls within the Si 
enriched layers for at least this period. By day +18 the blue wing has shifted significantly to 
the red, while the red wing remains constant except for the effects of an Fe II blend around 
6250 A. Other Fe II features begin to appear or strengthen at this phase as well. This 
behavior could be a signal that the Si layers are becoming detached from the photosphere 
by day +18. 

The day +25 spectra show a double minimum at the location of the Si II A6355 feature 
(see figure [7]). Telluric absorption is weak in this wavelength range, and the line profile is 
clearly seen in each of the three individual exposures, which support the reality of this feature. 
A possible explanation for this feature is contamination from the host that is not removed 
by the template subtraction; however, galaxy spectra do not typically exhibit features in this 
range that could cause such interference, and even if such were the case, we would expect 
to see similar behavior in the +34 day spectra. A seco nd possibility is contarnination from 



Fe II hues. Using the spectral analysis tool SYNOW (jjeffery &: Branch 



1990 



Fisher et al. 



1997 . 1999), and the example of SN 1994D as a starting place ( Branch et al. 20051 ). we find 



that while Fe II likely produces the absorption dips ~ 100 A away on either side of the Si II 
A6355 line, it is unlikely responsible for the double minimum. The third possibility, which we 
favor, is that this double minimum simply appears because we are resolving the Si II A6355 
doublet. This result implies that the Si II seen in the +25 day spectra is confined to a very 
narrow region of velocity space {Av ^ 1, 500 km s~^). If accurate, the true minimum of the 
Si II A6355 doublet would be about 100-200 km s~^ faster than indicated in figure [6] and table 
|2l but still significantly below the plateau velocity. The emergence of this thin layer may 
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also be responsible for the appearance of the narrow core in the +10 day spectrum as well as 
the apparent double minimum to the +18 day data. Some remnant of the blue component 
to the doublet may persist to the +34 day spectrum. Figure [7] also shows the spectra of 
several other SNe la taken around 25 days after maximum light. While the distinctly double 
minimum appears unique to SN 2005hj, the width and depth of the Si II feature is roughly 
consistent with the others. 

SN 20 05hi clearly b e longs to the low velocity gradient (LVG) group in the classification 
scheme of iBenetti et al.l ( 120051 ). but moreover the velocity derivative from maximum light 
through day +18, v = 3±7 km s~^ day~^, is consistent with no chang^. From the line profile 
evolution (table [21 figures [5] and [6]) we can deduce a plateau phase starting at —2.5 ± 2.5 
days which lasts no more than 30 days. Noting the change in the Si II A6355 line profile 
in the +18 day spectrum, we conservatively mark the end of the plateau phase as day 
17.5 ± 7.5 days, which gives the plateau phase a total duration of 20 ± 10 days. The Si II 
A6355 velocity evolution derived from the minima of FT smoothed spectra of several selected 
SNe la is plotted in figure [81 The velocity plateau of SN 2005hi is similar to that of other 



over-lu minous SNe la such as SN 1999aa (iGaravini et al.l l2004l ) an d SN 2000cx (Li et al. 



2001al ). bu t it is distinct from normal SNe la such as SN 1994D (jPatat et al.l Il996l ) and 
SN 1992 A ( Kirshner et al. 1993 ) that do not show a plateau phas^. 



3. Physical Constraints from Explosion Models 



There is general agreement that S Ne la result from some process involving the combus- 
tion of a degenerate C/0 white dwarf (iHoyle fc Fowled 1 19601 ). Within this general picture, 
two classes of models are most likely. The first is an explosion of a C/0- WD with a mass close 
to the Chandrasekhar limit {Mch) that accretes m atter through Roche-lobe overfiow from 
an evolved companion star (IWhelan &: Ibenlll973l ). In this case, the explosion is triggered 
by compressional heating near the WD center. Alternatively, the SN could be an explosion 
of a rotating configu ration formed from the merging of two low-mass WDs, after the loss of 
angular momentum (llben fc Tutukovl Il984j : [Paczynskil Il985l ) . Candidate progenitor systems 
have been observed for both scenarios: WD binary systems with the correct period to merge 



■^As defined bv iBenetti et al.l (|2005r ). v is the average daily rate of decrease in the expansion velocity 
from maximum light through the last available spectrum before the Si II A6355 line disappears; therefore, 
including the day +25 spectrum, SN 2005hj formally has v 
of freedom of 3.2 



27 + 4 km s ^ day ^ , but with a per degree 



''^We removed the data point at day —3 from the SN 1999aa curve because all spectral features in these 
data seem to be frequency shifted including the telluric features. 
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in an appropriate ti me scale with an appropriate total mass (IMaxted et al.ll2000l'): and super- 



soft X-ray sources (IGreiner et al.lll991 



van den Heuvel et al.l 



1992 



Rappaport et al.lll994 



Kahabka &: van den Heuvellll997l ) showing accretion onto the WD from an evolved compan- 
ion. There ar e still open questions about the details of both the merging and accret ion 
processes (e.g. lNomotdll982l : iBenz et al.lll990l : IPiersanti et al.ll2003l : iNomoto et al.l 120031 ). 



From the observed spectral and light curve properties, the first scenario appears to be 
the most likely candi date for the ma j ority of normal SNe l a. In particular, delayed deto- 



nation (DD) models (IKhokhlovl Il99ll : lYamaoka et al.lll992l : IWooslev fc Weaverl Il994j ) have 



and spectra of SNe la reasonably well 


(Hofiich 


1995; 


Fisher et a 


1995; 


Hofiich & Khokhlov 


1996; 


Lentz et al. 


2001; 


Hoflich et al. 


2002; 


Marion et al. 


2003 


1 i 


2006). 


In the DD scenario. 



a slow deflagration front turns into a detonation. The WD pre-expands during the deflagra- 
tion and undergoes com plete burning during the detonation phase. Similarly, the classical 
deflagration models W7 (INomoto et al.lll984j ) show similar behavior to DDs but only by ne- 
gle cting ins t abilit ie s due to the deflagration fronts (IGamezo et al.ll2003l ). For recent reviews 



see 



BranchI fll998h : iHillebrandt fc Niemeveil fl2000h : iHoflichI fl2006l ) 



Despite the success of classical DD and W7 models, both lack the basic features seen 
in SN 2005hj. Neither predicts a long plateau in velocity; they instead show a smooth 
decline of the photospheric velocitj]^ as a function of time (Fig. [9]). This happens because 
in expanding envelopes the photosphere recedes in mass and, because of the homologous 
expansion, in velocity as well. This behavior results from the smoothly declining density 
structure of the WD and the fact that variations in the speciflc energy production are small. 

In contrast, shell- like density structures will produce velocity plateaus i n a natural way 



becau se the photosphere remains in the shell for some time as shown by iKhokhlov et al. 



( 119931 ): iHoflich fc Khokhlovl (119961 ). To form a shell-like structure requires interaction of 
rapidly expanding material with a surrounding envelope. Various mechanisms have been 
suggested to supply this surrounding matter: the pulsating delayed detonation scenario 
( IHoflich et al.lll996l ). mergers or tamped detonation models. Shells n aay also form by the 



interaction of an ex ploding single WD within the progenitor system ( IGerardy et al.l 12004 
Quimbv et aDl2006h . 



^Although the Si II A6355 line is an imperfect tracer of the photospheric velocity as mentioned in i |2.11 the 
observed SNe la population typically exhibits a 1000-300 km s~^ decrease i n the measured line velocities 
between 1 week before maximum light to 2 weeks after ( Branch et al. 1988 : Benetti et al. 2005h . and the 
deflagration and classical delayed detonation models employed to explain t hese events have show n a cor- 



respondingly large decrease in photospheric velocities over the same period ([Khokhlov et al 
models are inconsistent with the corresponding w 300 km s^^ shift measured for SN 2005hj. 



1993h . These 
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We analyzed the observations of SN 2005hj based on detailed, spherical models for 
supernovae published in the hterature. The models are based on detailed calculations for 
the explosion, light curve and spectra. The models considered include delayed detonations, 
deflagrations, pulsating delayed detonations and tamped detonation/merger scenarios. In 
figure [9], we show the photospheric velocities as a function of time for these models along with 
the Branch-normal SNe la to illustrate the formation of a plateau in the models that naturally 
form a shell. Note for lower shell masses, this "plateau" is more accurately described as a 
period of slowly declining velocities. 

In classical delayed detonation models and for normal-bright SNe la. Si is present over 
a wide range in mass, spanning about 0.4 to O.SM©, which corresponds to velocities from 
about 8,000 to 9,000 km s~^ to more than 20,000 km s~^. The Si layer is thick (in the mass 
frame) because explosive oxygen burning occurs over a wide range of temperatures. The 
density gradient is smooth and Si is mostly in Si II, so initially the velocity derived from the 
minimum of the Si II A6355 line smoothly declines with the receding photosphere governed 
by the geometrical dilution of the expanding envelope. Eventually, the photosphere begins 
to recede below the Si layer at which point the evolution of the Si II A6355 line profile will 
show the following behavior: 1) the optical depth of the highest velocity material at the 
largest radii will begin to decline below 1 and as such the blue wing of the line profile will 
start to drift toward the red; 2) as the optical depth decreases, the strength of the line as 
measured from the line depth will decrease; 3) the line minimum may continue to slow, but it 
will grow increasingly discrepant with the photospheric velocity. This phase typically begins 
1-2 weeks after maximum light for normal SNe la and is heralded by the appearance of Fe II 
lines. While this behavior is commensurate with observations of normal-bright SNe la such 
as SN 1994D, this behavior is not consistent with the observations of SN 2005hj. 

The Si II A6355 line seen in SN 2005hj is narrow, and during the plateau phase the 
wings do not change, the depth does not change, and the velocity of the minimum does not 
change to within the errors. The data require a narrow distribution of Si II in velocity space, 
and we suggest this may be explained by an interaction that compresses the Si rich layers 
as predicted by merger and pulsating delayed detonation models. The shell models are also 
consistent with the velo city drop seen after the plateau because a signific ant amount of Si is 



located below the shell flKhokhlov et al.lll993l : iHoflich fc Khokhlovlll996[ ) 



In Fig. [To], we show general properties of these models. As discussed in the papers 
above, to first order, the observational signatures of the shell depend on the mass of the 
shell M shell- Almost the entire WD is burned, and momentum conservation determines the 
amount of high velocity matter that can pile up in the expanding shell. With increasing shell 
mass, more material of the SN envelope is slowed down. As a consequence, the velocity of 
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the shell Vsheii decreases with M shell- Because it will take longer for the photosphere to recede 
into the lower velocity matter, the time until the beginning of the plateau phase, to? increases 
with M shell- The optical depth increases with Msheiu duration of the plateau. At shell, also 
increases, the temperature gradient becomes st eeper, and the photosphere becomes cooler 
(i.e. B — V increases) with increasing Msheii (IKhokhlov et al.l Il993l : iHoflich &: Khokhlov 
19961 ). The duration of the plateau, Atsheiu is defined by the velocity spread 5v around 
V Shell with 5v = 500 km s~^, which puts the end of the plateau phase safely into the parts 
of a rapidly declining Vph- We choose a larger value than in the observations to avoid 
ambivalences due to discreetness, which, in some of the models, is of the order of ~ 100 km 
s"^. By increasing 6v from 200 to 500 km s~^ the nominal duration is increased by ^ 1 day. 
However, we also note that the actual width depends on the velocity spread in the shell (see 



Given the model predictions, we can use different observational indicators to test which 
M shell is consistent with SN 2005hj (Fig. [101). three parameters, Vsheiu t shell, and to 
suggest Msheii ~ O.2M0, with the allowed ranges specifically bracketed by 0.15-0.6, 0.1-0.25, 
and 0.1-0.25 M© for the plateau length, shell velocity, and plateau onset, respectively, taking 
the observed errors into account. The comparison between the B — V color as a function 
of V shell, i shell, or Iq, however, shows only marginal consistency between the observations 
and the models if we assume only foreground redding by the Galaxy. We note, that the 
intrinsic B — V color of the models is uncertain by about 0.05 to 0.1™ at maximum light. 
The two best-fitting models, pdd3 and det2env2, show a peak brightness, My, of —1 9.42'" 
and —19.41"*, respectively, with an uncertainty of 0.1"* (IHoflich fc Khokhlovlll996l ) vs. a 
typical DD model with -19.2™ jHofiich et al.ll2002h . i.e. they are brighter by about 20 % 
mostly due lower escape probabl ity of 7-rays that results when the ^^Ni layers are slowed 
down because of the interaction (jHofiich et al.lll99ll ). 



4. Discussion and Conclusions 

We have presented photometric and spectroscopic data for SN 2005hj, a slightly over- 
luminous Type la. The most striking feature is an apparent plateau in the expansion velocity 
evolution, which we derive from the Si II A6355 line. The velocities remain at about 10, 600 
km s~^ for about 3 weeks starting slightly before maximum light, and this plateau is bracketed 
by preceding and succeeding decelerations. We find that Si is confined to a relatively narrow 
velocity region. Analysis of the detailed observations in concert with published models 
suggest there may be some physical distinction between SN 2005hj and other normal-bright 
SNe la that may systematically affect their use as distance indicators if not properly taken 
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into account. 



The models considered include delayed detonations, deflagrations, pulsating delayed 
detonations and tamped detonation/merger scenarios. In order to explain the narrow Si II 
A6355 line and its plateau in velocity, we suggest an early interaction that forms a dense shell 
as predicted by merger and PDD models. The spectral and photometric peculiarities are 
consistent with respect to the velocity, duration, and onset of the plateau, and marginally 
consistent with the maximu m light color, for models that have shells of abo ut 0.2Mq. As 
indicated by earlier works (IKhokhlov et al.l Il993l : iHoflich fc Khokhlovl 119961 ) , the mass of 
the interacting shell has been found to be the parameter that dominates the details of these 
observational signatures independent of how this shell may form. The tight predicted relation 
between each of v shell, t shell, and to provide a stable means to separate SN 2005hj like 
events from regular Branch-normal SNe la. Although the agreement between the shell models 
and the observations is good, the predictions are not necessarily unique and other possibihties 
may exist. For example, we have not considered 3 -D models such as the detonation from 
a failed deflagration scenario recently examined by iKasen fc Plewal (120061 ). For SN 2005hj 
then the agreement of the plateau velocity and its duration to that predicted by shell models 
may simply be a fluke, and in such case this concordance should then not hold for other SNe 
with similar Si II A6355 evolution. Given the data and models considered, we suggest either 
PDDs or merger events are responsible for SN 2005hj, and this implies the existence of two 
different progenitor channels. 

It is important to understand how these two progenitor channels, which may occur in 
relatively varying f ractions a s a fun ction of redshift, will impact studies using SNe la as 
distance indicators. iLi et al.l (j2001bl ) estimate that 20% of SNe la in their sample are either 
1991T-like or 1999aa-l i ke. T hese SNe show spectral features a nd a velocity pla teau similar 
to SN 2005hj. iBranchI J200l[ l found 5 1999aa-like events in the lli et all J2001b[ l sample out 
of 20 total SNe la that were observed early enough to show 1999aa-hke spectral features, 
and one that was 1991T-like; however, in the pre-LOSS sample they do not classify any 
of the 7 SNe la with early spectra as 1999aa-like. These nearby samples are constructed 
from targeted galaxy searches that have different selection biases than the deep cosmological 
surveys, but we will assume a uniform 2005hj-like rate of 25% for all SNe la. SNe la that 
appear spectroscopically similar to SN 2005hj in a single epoch could none the less arise from 
different progenitors, and the mass of the low-density envelope around PDDs or mergers may 
effect their peak magnitudes and/or light curve shapes, but we will fur ther assume that al l 
such events deviate uniformly from the LWR of Branch normal SNe la. iHoflich et al.l (119961 ) 
calculated the relation between peak V band magnitudes. My, and the fading between 
maximum light and -|-20 days, dMy (20), for a variety of theoretical models and found that 
shell models produced dMv{20) that were 0.2 to 0.3 mag smaller than for (in vacuum) 
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delayed detonations reaching the same peak magnitude. Therefore using the same LWR for 
shell models will result in corrected peak magnitudes systematically offset by 0.1 to 0.2 mag. 
Also, the observed peak magnitudes of SNe la are usually corrected for absorption along the 
line of sight using the observed B — V color at maximum light and a reddening law. For 
events that are intrinsically red, this will increase the estimated peak magnitude above its 
already over- luminous intrinsic value. Cosmological studies may therefore need to remove or 
at least separately handle SN 2005hj-like events to avoid systematic errors in distance. 

As a case for the importance of separating different progenitors, let us consider SN 1999ee. 
Very similar to SN 2005hj, SN 1999ee shows a plateau with Vsheii = 1 0, 500 km sec~\ a d ura- 
tion of 14 ±3 days, and an onset at day — 3± 1 relative to maximum (jHamuy et al. 2002 : see 
Fig. Ej- The B-V color of SN 1999e e was also quite red at maximum light; +0.28 "" ± 0.04 



after correction for galactic extinction (jStritzinger et al.ll2002l : iKrisciunas et al. 



20041). Based 



on the standard brightness decline relation and the corresponding colors, iHamuy et al.l (120021 ) 
derived reddening in the host galaxy of 0.28 ± .04 which implies an absolute brightness of 
Mv = -19.95™ similar to SNLS-03D3bb which iHowell et al.l J2006h attributed to a super- 
Chandrasekhar mass WD. Taking into account the spectroscopic information about the 
velocity plateau, its length and onset, we attribute a portion of the red color to the intrinsic 
properties of the supernova. We find that the duration of the velocity plateau, its onset and 
size are consistent with a shell mass of O.2M0 which suggests an intrinsic color B — V of 
0.15™ ± 0.02 (see Fig. [TOl) . This reduces the reddening in the host galaxy to ~ 0.13™ and 
the absolute brightness My to —19.53™, which compares favorably to the model predictions 
of —19.42™ and —19.41™ for pdd3 and detenv2, respectively, within the model uncertain- 
ties. Note that there is an interstellar sodium line in the spectra that implies some degree 
of reddening within the host. There are some apparent spectral differences when compared 
to SN 2005hj, namely SN 1999ee has a slightly broader blue wing in Si II and stronger 
absorption around 4900 A. This may either hint toward different explosions scenarios (i.e. 
pulsations versus mergers), or different viewing angles of asymmetric envelopes. 

This brings us to the limitation of our studies. Except for the color, SN 2005hj fits 
remarkably well with the merger and FDD model predictions but, still, it is a single event 
and the good agreement may be coincidental. We need a large, well-observed sample of 
similar objects to test and verify or falsify the models and to determine the shell mass 
distribution. Moreover, 3D effects have been neglected. In reality, we must expect some 
dispersion. Though pulsating delayed detonation models may be expected to be rather 
spherical, mergers may be asymmetric with directionally dependent luminosities and colors. 
In fact, both classes may be realized in nature. As mentioned above, the duration of the 
plateau. At shell, is defined by the velocity spread around vsheii- The physical width of the 
shell depends, to first order, on the distance at which the interaction occurs and the density 
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distribution of the interacting expandin g media and shell during the hydrodynamical phase 
of the interaction ( Gerardy et al.l 120041 ) . For obvious reasons, asymmetries of the shell will 
increase the velocity gradient seen over the shell. The observations of SN 2005hj indicate 
a very flat plateau that, in principle, may further constrain the properties of the shell. 
For SN 2005hj, this may already indicate a rather spherical shell and hint toward the FDD 
scenario or mergers with an intermediate disk of very large scale heights. However, additional 
information needs to be taken into account such as detailed line profiles and statistical 
properties to break this degeneracy between mergers and PDDs. 

As a next step, detailed models for the explosion, light curves and spectra tailored 
toward individual events need to be constructed. Whereas the mean velocity of the shell for 
a given mass is dictated by momentum conservation, the thickness of the shell is limited by 
the distance of the shell material, the distance sound can travel during the interaction, and 
the specific density profile within the shell. With increasing distance of the shell, the relative 
size (and corresponding velocity spread) becomes smaller because the sound speed remains 
about the same. The intrinsic color will be sensitive to the optical depth of the shell, which 
is governed by the magnitude o f the density iump a nd thus depends on the distance of the 
interacting shell from the WD (IGerardy et al.ll2004j ). The blue B — V color for SN 2005hj 
may hint of a need to modify the distance and structure of the shell. Frecise analysis of such 
"non-stable" features requires detailed model fitting beyond the scope of this paper. 

In the recent past, both the scenarios leading to shell-like structures have been dis- 
counted. FDD models have been dismissed because 3D defl agration models showed that the 
WD becomes unbound and thus pulsations would not occur ( Ropke et al.ll2003 : Hillebrandt fc Niemeyer 
2000). However, it has recently been shown that this solution depends mainly on the ignition 
conditions, namely the number and locations of ignition points leading to single or multiple 
bubble solutions, and mixture of bubble solutions leading to Raleigh- Taylor instabilities. As 
a result, solutions with fewer bubbles are likely to result in a reduced amount of burning. 



thus only slightly unbin ding the WD and increasing the possibility of FDDs (ILivne et al. 



20051 : iFlewa et al.l 120041 ). Similarly, the merging scenario has been dismissed because the 



WD may undergo burni ng during the merger and result in an accretion induced collapse 
( iNomoto &: Kondd Il99ll ). and also on the basis of the long merging time scale. However 
both of these results depend sensitivel y on the initia l conditions, and new pathways to the 
actual merging may effect the results (iLii et al.ll2006l ). In light of our results, the predicted 
death of both of these scenarios may be premature, and further studies are needed. 
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Table 1. Observing Log for HET Spectra of SN 2005hj 



Date (UT) 


JD-2400000.5 


Phase (day) 


Exp. (s) 


Filter 


Slit (") 


Oct 26.30 


53669.30 


-6 


2x600 


GG385 


2.0 


Oct 27.20 


53670.20 


-5 


4x600 


OG590 


1.5 


Oct 27.28 


53670.28 


-5 


2x600 


GG385 


1.5 


Nov 1.27 


53675.27 





2x600 


GG385 


1.5 


Nov 1.29 


53675.29 





2x550 


OG515 


1.5 


Nov 3.27 


53677.27 


2 


2x480 


GG385 


1.5 


Nov 4.26 


53678.26 


3 


2x600 


OG515 


1.5 


Nov 6.17 


5.3680.17 


5 


2x600 


GG385 


1.5 


Nov 6.18 


5.3680.18 


5 


2x600 


OG515 


1.5 


Nov 11.25 


53685.25 


10 


3x600 


OG515 


1.5 


Nov 19.22 


53693.22 


18 


4x600 


OG515 


1.5 


Nov 26.20 


53700.20 


25 


2x600 


OG515 


1.5 


Dec 5.18 


53709.18 


34 


3x600 


OG515 


2.0 



Note. — Phases are relative to the derived date of maximum light, 2005 Nov. 
1.6, and rounded to the nearest day. Effective wavelength ranges for the GG385, 
OG515, and OG590 setups are 4100-8900 A, 5100-9800 A, and 5900-10,800 A, 
respectively. All spectra from a given night were combined for the figures and 
analysis. 
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Table 2. Derived Quantities for the SN 2005hj Si II A6355 Line 



Phase 
(days) 


Velocity 
(km s-l) 


""vol 

(km s^-'-) 


Depth" 


FWHM 

(A) 


-6 


10820 


140 


0.38 


140 


-5 


10800 


110 


0.32 


120 





10640 


90 


0.52 


110 


2 


10440 


100 


0.60 


110 


3 


10640 


90 


0.57 


110 


5 


10680 


80 


0.57 


110 


10 


10530 


100 


0.60 


100 


18 


10550 


120 


0.48 


90 


25 


9850 


90 


0.25 


60 


34 











Note. — Quantities measured from the FT smoothed 
data. The Si II A6355 line is not clearly detected in the 
+34 day spectrum and thus no measurements are reported. 

"Similar to [Leonard et al. I ||2002| '). depth is defined as 
(/c — /min)//c where /min is the flux at the minimum of 
the smoothed line and fc is the estimated continuum level 
at the corresponding wavelength. 
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17.01 I 




-20 -10 10 20 30 40 

Days After 2005 Nov 1.6 



Fig. 1.— ROT SE-IIIb unfil t ered l ight curve of SN 2005hj (filled circles). The best fit /?-band 
template from iKnop et al.l (120031 ) is plotted as a solid line over the fitting range, and as a 
dotted line continuing on to later phases when the rapid decline of the flux below ~5500 A 
causes our unfiltered light curve to fade faster than the i?-band decline rate. Arrows mark 
5-cr upper limits of the subtractions determined from the noise level in annuli centered on 
the location of the SN. Epochs with HET spectra are marked with "S." 
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Fig. 2. — HET/LRS spectrum of SN 2005hj at —6 days (black line). The blue curve shows 
the relative contribution from a SN 1999aa-like superno va at the same phase from a linear 
interpolation of the —7 and —3 day spectra presented by lGaravini et al.l (120041 ). and the red 
line gives the estimated relative host galaxy light contamination. The relative scaling for the 
galaxy and SN components were determined through a least squares fit to the SN 2005hj 
spectrum. The purple line shows the fit (portions of the spectra effected by telluric lines 
were excluded from the fitting as shown by gaps in the fit curve). 
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Fig. 3. — Spectral evolution of SN 20051ij recorded by the HET from —6 to +34 days after 
maximum light (2005 Nov. 1.6). The estimated galaxy contamination has been subtracted 
and the spectra have been shifted for clarity. For display purposes, portions of the spectra 
with low signal to noise have been smoothed (thin line segments). The typical telluric 
absorption spectrum is shown by the grey shading along the top of the figure. 
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Fig. 4. — Compar ison of the spectra of SN 2005hj (colored lines) and SN 1994D (black lines; 



Patat et al.lll996l ) at 5 days before and 2 days after maximum light. The spectra have been 
shifted for clarity. At both phases shown, the absorption component of the Si II A6355 line is 
more narrow for SN 2005hj. The relative line depths for this feature are similar at +2 days, 
however, the Si II A6355 absorption is much weaker for SN 2005hj in the —5 day spectrum. 
Despite these differences, the velocities inferred from the minima of the Si II A6355 lines are 
consistent between SN 2005hj and SN 1994D to within 5%, the former being 560 km s~^ 
slower at —5 days and 150 km s~^ slower at +2 days. 
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Fig. 5. — Evolution of the Si 11 A6355 line. In the top plot, the spectra have been shifted 
by arbitrary factors for clarity with the youngest spectrum on the top and the oldest on the 
bottom, and they are colored as in figure [3l Each spectral bin is plotted as a rectangle with 
the height indicating the statistical 1-a error interval. The effects of the FT smoothing are 
shown by the grey curves, and open circles mark the FT smoothed minima. The bottom 
plot shows the data normalized by the continua to emphasize the relative spectral evolution. 
The vertical dashed line marks the plateau velocity of 10,600 km s~^. 
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Fig. 6. — Velocity of the Si II A6355 line minimum as determined from the FT smoothed 
minima. The error bars are the quadrature sum of the 1-a confidence intervals found via 
Monte Carlo simulations and the 80 km s~^ scatter found in the H II region line centers 
among the different epochs. 
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Fig. 7. — Spectra near the Si II A6 355 doublet for S N 2005hi and other SNe la ar ound 25 
days after maximu m light f l994D: IPatat et al.lllQod : 1999aa: [Garavini et allbooi 1999ee; 
Hamuy et al.ll2002l : 2003cg; lElias-Rosa et al.l l2006l ) . Vertical dotted lines mark the doublet 
components blue shifted by 10,000 km s~^. The doublet is clearly resolved only in the case 
of SN 2005hj, while the fiat bottom to the SN 1999aa profile may suggest similar behavior. 
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Veloci ty of the Si II A6355 line mi nima for SN 20051ij and oth er well-observed 



SNe la ri992A: iKh^shner et all Il993l: 1994D 



2000cx: 



Li et al 



200 lal : 1999aa; 



Patat et al. 1996 



Garavini et alj 



1998bu; 



Jha et al. 1999 



20041 : 1999ee; iHamuv et all l2002l : 2005cg; 



' ^ 'I I ' ■ ' 'I I ■ yi ' ■ ' 'I <~> I 

Quimby et al.ll2006l ). All velocities were calculated from the FT smoothed minima using the 



same filter. The ordinate shows the change in velocity relative to the (interpolated) day +20 
value. Lines connecting the data points are intended only to guide the eye. The sample was 
not drawn randomly from the parent population; instead it was tailored to show the range 
of behaviors present in well-observed in normal-bright events. 
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Fig. 9. — Photospheric velocity versus time since explosion for models of Branch-normal 
SNe la. The photosphere is defined by the radius of last scattering, and because Thomson 
scattering is the dominate source of continuum opacity, this position is largel y wavelength 



indep endent. Shown are the classical delayed detonation model 5p 0z02.25 (IHoflich et al. 



2002h. the pulsatin g delayed detonation models pdd3 and pdd9 flKhokhlov et al 



1993 



Hoflich et al.l Il995bl ). and the tamped detonation models det2en2/4 ( iHoflich fc Khokhlov 



19961 ). The X's mark the velocities of SN 2005hj as measured from the minima of the Si II 



A6355 line assuming a 22 day rise to peak. The dotted horizontal lines mark the region with 
> 1% Si abundance in model pdd3. As discussed in the text, the Si II A6355 line minimum 
is an imperfect indicator of the photospheric velocities, and a systematic offset can develop 
at later times once the photosphere has receded below the Si rich layers. This may explain 
the discrepancy between the model predictions and the last two observations, which show 
weakening Si and strengthening Fe lines. The low velocities observed at early times for 
SN 2005hj may indicate that rather little outward mixing took place. 
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Fig. 10. — Properties of Branch-normal SNe la models with interacting shells of mass Msheii 
including delayed d etonations (5p0z22.25 with Msheii = 0.01,0.05 and O.l Mf;^, open circles 
Gerardv et al.lbooi). pulsating d elayed detonation raodels ( open stars; pdd3; iKhokhlov et al.l 



1993; 



Hoflich fc Khokhlovlll996l: pdd9: iHofiich et al.lll995br) and merge r models (open trian- 



gles; det2env2/4/6; iKhokhlov et al.lll993l : iHoflich fc Khokhlovl Il996l ). In the upper plots, 
we give as a function of Msheii the duration of the velocity plateau Atpiateau, its velocity 
Vpiateau defined by 5v < 500 km s~^ (see text), and its starting point to relative to maximum 
light. In the lower plots we give observational relations, namely B — V aX maximum light as 
a function of Vpiateau, ^tpiateau and to- Note that B — V is uncertain by about 0.05 to 0.1™. 
In addition, we give the observations of SN 2005hj (filled squares) including error bars. 



